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THE COMPRESSIBLE LAMINAR BOUNDARY LAYER WITH HEAT TRANSFER AND ARBITRARY 

PRESSURE GRADIENT ' 

' By Clabence B. Cohen and Eli Resiiotko 


SUMMARY 

- An approximate method for the calculation of the com- 
pressible laminar boundary layer with heat transfer and 
arbitrary pressure yradient, based on Thwaites’ correlation 
concept, is presented. The method results from the application 
of Stewartson’s transformation to Prandtl’s equations, which 
yields a nonlinear set of two first-order differential equations. 
These equations are then expressed in terms of dimensionless 
parameters related to the wall shear, tjie surface heat transfer, 
and the transformed free-stream velocity. Thwaites' concept 
of the unique interdependence of these parameters is assumed. 
^The evaluation of these quantities is then carried out by utilizing 
exact solutions recently obtained. 

With the resulting relations, methods are derived for the 
calculation of the two-dimensional and axially symmetric 
laminar boundary layer with arbitrary free-stream velocity 
distribution, Mach number, and surface temperature level. 

The combined effect of heat transfer and pressure gradient 
is demonstrated by applying the method to calculate the char- 
acteristics of the boundary layer on thin supersonic surfaces 
and in a highly cooled, convergent-divergent, axially symmetric 
rocket nozzle. 

INTRODUCTION 

j In recent years, with the advent of laminar airfoils and 
with the observation of laminar boundary layers at Keynolds 
numbers as high as 50X10® (ref. 2), the ability to reliably 
estimate viscous flow and heat-transfer effects for a laminar 
boundary layer has become increasingly important. More- 
over, with high-altitude flight becoming more common, the 
subsequent lower Ecynolds numbers encountered should 
more frequently produce a laminar boundary layer. Sta- 
bility calculations based on the theory of Lees and Lin 
(ref. 3) have also emphasized the possibilities of maintaining 
a laminar boundary layer through cooling of aerodynamic 
surfaces. The effect of favorable pressure gradients in 
increasing the stability of laminar boundary layers may 
also make solutions to the laminar problem applicable to 
the design of nozzles and turbine blades. 

’ Solutions of the laminar-boundary-laycr equations that 
include effects of compressibility, pressure gradient, and 
beat transfer have been quite limited in number. Of the 
exact solutions, most have restrictions of range or applica- 
tion, or both. The solutions of references 4 and 5 are re- 


stricted to zero pressure gradient, while those of reference 0 
allow small pressure gradients. The developments of refer- 
ence 7 are restricted to small beat transfer and low Mach 
number. Solutions obtained by assuming that fluid proper- 
ties are constant or that the Mach number is essentially zero 
are obtained in references 8 to 10. Those solutions of 
references 11 to 13 that arc for a Prandtl number of 1 are 
not restricted in range of compressibility, pressure gradient, 
or beat transfer. However, they apply to specific types of 
free-stream velocity distribution that are inappropriate for 
general practical problems. 

In 1921, von Ktirman (ref. 14) recognized that to solve the 
skin-friction problem it was not necessary to have the exact 
and complicated solution, but that it would be quite satis- 
factory to evaluate average quantities across the layer if 
they could bo related to the surface values. The concepts 
of displacement and momentum thicknesses were introduced, 
thus considerably simplifying the mathematics of the prob- 
lem. With this integral method, if the form of the velocity 
profile is related to a single parameter, a method of calculat- 
ing the boundary layer is obtained. Pohlhausen (ref. 15) 
carried out this method by postulating a quartic velocity 
profile depending upon the local pressure gradient. A num- 
ber of investigators have extended Poblhau^n’s method to 
compressible flows over insulated surfaces. 

With the presence of heat transfer at the surface, the com- 
pressible problem becomes more complex. Kalikhman (ref. 

, 16) defined certain heat-flow quantities analogous to the 
^displacement aiid momentum thicknesses and, in a manner 
^ similar to Pohlhauscn’s, developed a complex iterative 
‘ procedure for the solution of the general problem. ' More 
' recently, references 17 to 20 have further developed this 
technique. The preceding methods arc tedious, since they 
require a solution of at least one ordinary differential equa- 
tion for any particular problem. 

Thwaites’ method (ref. 21) docs not requii-o the solution of 
ordinary differential equations. In that formulation, it is 
suggested that the basic goal of an integral approach might 
he achieved if the problem is considered as that of relating 
the wall shear, its normal derivative at the wall, and the 
form factor (ratio of displacement thickness) to one another 
withoutspecifyingatypeofprofile. To this end, nondimension- 
al forms of these quantities were defined and were evaluated 


I 1 Supsrsedes NACA TN 3325, “The Compressible Laminar Boundary Layer with Heat Transfer and Arbitrary Pressure Gradient,” by Clarence B. Cohen and Ell Re hotko, 1955. 
i 2 The principal developments of this paper, which Is part of the doctoral dissertation of the senior author (ref. 1), were carried out under the stimulus and guidance of Pro essor Lu g 
drocco and the sponsorship of the Daniel and Florence Guggenheim Foundation. The final analysis and the computations were completed at the NACA Lewis laboratory. 
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by examining exact solutions for the incomi)ressible laminar 
boundary layer. It developed that a nearly' universal 
relation existed between these quantities for favorable 
pressure gradients, and for adverse pressure gradients 
Thwaites selected a single representative relation. A unique 
correlation was chosen that reduced the solution of an incom- 
pressible problem to the evaluation of a single integral.® 

Rott and Crabtree (ref. 23) recognized that, in the absence 
of heat transfer and with a Prandtl number of 1, the Illing- 
worth-Stewartson correlation between compressible and 
incompressible boundary-layer solutions (ref. 24) could be 
used to extend Thwaites’ method to include effects of 
compressibility. 

With the presence of heat transfer the application of 
Stewartson’s transformation docs not correlate a given com- 
pressible problem to an equivalent incompressible one. 
Thus, the universal relation previously described is not 
adequate. Unfortunately, there is little possibility of 
establishing a family of universal relations with, for examiilc, 
the wall temperature as the distinguishing parameter, since a 
variety of exact solutions to this problem is not available. 
However, one such sot of relations may bo obtained from the 
solutions of references 11 to 13. 

In the present paper, after formulation of a nonlinear 
system of two first-order differential equations (with the 
major restriction being a linear viscosity law), methods of 
solution arc developed depending on Thwaites’ concept of 
universal functions. The functions used for this purpose 
are evaluated from the solutions of reference 12 only. 




BOUNDARY-LAYER EQUATIONS 


The equations of the steady, two-dimensional compressible 
laminar boundary layer for perfect fluids are 
Continuity: 


Momentum: 


^(p«)+|^(pi’)=0 


du , du bp , d 
bx ^ by~ bx b7j 




Energy: 


by 


=0 


( 1 ) 


( 2 ) 


dA . bh bp , b / fi bh\ I /duV 
^ bx'^^'^’ by~'^ bx'^by \Pr by)'^'" \by) ® 


flow. The constant X is used to match the viscosity with 
tlic Sutherland value at a desired location. If this location 
is assumed to be the surface, the result is 



where 


^,„= Sutherland’s constant= 198.0° R (for air) 
Stewartson’s transformation. — The velocities in the equa- 
tions of motion (1) to (3) can be replaced through tlie 
definition of a stream function: 





Introducing the quantity X from equation (4), a slight modi- 
fication of Stewartson’s transformation (ref. 24) may bo 
written: 



( 7 ) 


The transformed coordinates are now represented by upper- 
case letters {X, Y), and the subscript e refers to local con- 
ditions at the outer edge of the boundary layer (external). 
The subscript 0 refers to free-stream stagnation values. 

Applying equations (4) to (7) to the boundary-layer 
equations (1) to (3) and assuming that Pr and Cp are constant 
(but not yet requiring that Pr=l) result in the following 
system: 

f7x+I'r=0 • (8) 

UU:c+VVy=U,Ue^{l+S) + y„Ury (9) 



■ • ( 10 ) 
where the enthalpy term S is defined for convenience as 


(All symbols arc defined in appendix A.) 
The viscosity law to bo assumed is 



( 11 ) 


Equation (4) is of the form taken by Chapman and Rubesin 
(ref. 5), except that the reference conditions (yo, Q are free- 
strcarfi stagnation values, since, in' the case of pressure 
gradient, the local stream values are not constant along the 


where h, is the local stagnation enthalpy. The stream 
function has been replaced by the transformed velocities 
( U, V) through the relations 


U=i/y 


( 12 ) 


8 Other approaches, such as that of Young and 'VVinterbottom (ref. 22), have resulted In 
expressions for the momentum thickness similar to that of Thwaites. In that analysis, 
however, the derivation was a modification of the Pohlhausen technique. The application 
of a correlation concept was not proposed. 


The resulting relation between tbe transformed and physical 

longitudinal velocities is U=— u. 
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I The boundary conditions applicable to the system (8) to 
(10) for a specified wall temperature are 


f/(X,0)=0 

F(Z,0)-0 

SiX,0)=SM) 

lim)S'=0 


(13) 


Ihn U^^ileiX) 

Y—*<x> 


, Integral equations. — An alternate form of the momentum 
equation may bo obtained by subtraoting the momentum 
j equation (9) from the product of the continuity equation 
1(8) and the quantity (17^— U). This results in 


:[U{u-u)]^:^[V(u,-u)]y+u,^(u-U)+ 


UeU,^S=- 


■ro 


Uyy (14) 


■ If equation (14) is integrated with respect to Y between the 
i limits y=0 and Y=A, where A is a constant distance normal 
I to the surface sufficiently large that the conditions S=0 
(and U=Ue can both be satisfied, there results 


d 


^ ienUD + U,UeJ1r = MUy)y=, 


(15) 


(Where the transformed momentum thickness fl,, and the 
transformed displacement thickness are defined as 




(16) 


I By carrying out the indicated differentiation, equation (15) 
can be put in the form 


U, 


(20<,+3f,)= 


dX 


Ul 


{Uy\ 


(17) 


I This equation has the form of the conventional Kerman mo- 
' mentum integral. 

I It should be noted that because of Stewartson’s trans- 
1 formation a simple relation exists between the parameter 
■ dtr and the actual physical momentum thickness d. This 
relation is 

I 7+1 


\2(7-U 


(18) 


Following a procedure with the energy equation similar to 
that for the momentum equation results in 


dE U^E- 

dX^V, 


Vo 


PrU, 




Where the convection thickness is defined by 

U' 


-j: 


Sjj-dY 
^ 0 . ■ 


(19) 


(20) 


The method presented in this report uses exadt solutions 
to the boundary-layer equations including tfie energy equa- 
tion. Since both the skin-friction and heat-transfer param- 
eters from the exact solutions are correlated with a parameter 
which is evaluated from only the momentum intcgi-al equa- 
tion, it will not be necessary to use equations (19) and (20). 

REDUCED BOUNDARY-LAYER EQUATIONS 

At this stage, the relation governing the boundary-layer 
Icvelopmeut is equation (17), subject to the bouiidary condi- 
dons {d,r)T=o=0 or (g<,),.,. where the subscript sp indicate s 
stagnation-point valu es. The former condition on 6,r applies 
when the boundaryTayer is initiated without a stagnation 
point (such as in the case of a supersonic thin airfoil). The 
value of {6,r)gp depends on the value of (Uey)sp and on the 
surface temperature. At a stagnation point; {Btr)sv—^sp- 
'Vnlues+)fr 0 , 7 TircrprcscntoduTutablcnl . 

Before consideration of a solution that depeiids on a cor- 
relation similar to that of Thwaites, it is e.xpedient to trans- 
form the preceding system of equations to a system involving 
dimensionless parameters. The correlation concept will 
then be introduced and methods of solution developed. 

.||i(.^parametric form 

The dimensionless parameters, which arc related tojxsrriis^ 
appearing in equations (17) and (19), can be defined and 
evaluated from the following expressions: 

Shear parameter, or first-derivative parameter. 

Correlation number (related to pressure gradient), or second- 
derivative parameter, ( 


Vn 


/W 

eir 


). 


\te J \^e/ 


( 22 ) 


Ue 

Heat-transfer parameter, or third-derivative parameter, 

(&-).=”» (01. 

In definitions (22) and (23), use is made of the following 
relations, respectively: 


UelJ. 




(id-iSa) 


which is obtained from equation (9), and 
/^\ _■ 


(24) 


(25) 


which is obtained by differentiating equation (9) with respect 
to Y and evaluating the rcsidting expression at F=0. 


If equation (17) is multiplied by 


29>rUe 

Vo 


there rcsidts 


Ue d{e]r) 
Vo dX' 






■if®!! 


where Tit,—-— is the form factor for low-speed flow {Me 

v[T 


(26) 


= 0 ). 




I 
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TABLE I.— INITIAL VALUES OF PARAMETERS 
(a) Stagnation-point flow I 



C/-J1EZ 


Nuf’yflieZ 


fii7- /< 

A 








(7 



ir.r 


Pr=l 

Pr^O.7 

Pr = l 

Pr=0.7 

/V=0.7 

\ Vw ' * 


\ Vv, 

Prr^l 



(•) 


C>) 

(•) 

/V = l 


Pr=l 


Two-dimensional (^=1) 

- 1.0 

«1.30 * 

1.21 

c 0. 500 

0. 438 

0. 430 

0 -0.170 


«0. 400 

c -0.425 ; 

-.8 

« l.fiO 

1.49 

« . 522 

.452 

.451 

«■ .012 


® . 384 

O.031 

-.4 

0 2. 04 

2. 00 

® . 540 

.474 

.475 

.*345 


» .338 

« 1.021 

0 

2. 40 

2. 40 

.570 

. 495 

.495 

.648 


.292 

2.218 

1.0 

3. 47 

3. M 

.015 

.5:43 

.537 

1.386 • 


.177 

7. 850 

Axtiilly symmetric ($='A) 


- 1.0 

1.04 

1.50 

0. 700 = 

0.607 

0. 607 

-0. 0771 


0 . 3(H) 

-0. 257 

-.8 

1 . 8 f) 

1.79 

.712 

.617 

.621 

. 0576 


.289 

.199 

-.4 

2 . 2 ri * 

2. 22 

.739 

. 639 

. 643 

.318 


. 209 

1.185 

0 

2 . 02 

2. 62 

.76:1 

.662 

. 662 

. 569 


.248 

2. 298 

1.0 

3. 49 

3. 55 

.809 

.701 

.708 

1.105 


.194 

0.012 


* Thes?o values arc obtained when ens. (IR) of ref. 12 arc solved for Pr=0.7, JV/«=0. 
t These values arc obtained by multiplying the results for Pr^l by (0.7) 

« Interpolated values from solutions of ref. 12 . 

■ (b) Sharp edge or pointed body 


V 



Nu}-yj Re w 


Pr=l . 

Pr=0.7 

C^) 

Two-dimensional 

0. 664 

0.332 

0. 295 

Axially symmetric 

1. 150 

0. 575 

0.510 


• These values are obtained by multiplying 
the results for /^*»1 by ( 0 . 7 )>/ 3 , 

TABLE IL— SyAfMARY OF HEAT-TRANSFER AND WALL- 
SHEAR PARAMETERS | 


s„ 

0 

n 

1 

N 

T 

/CfRe^\ ' 

\ Nu. J e,-i 

Iltr 

- 1.0 

-0. 326 

0. 1335 

0 

1.0845 

0. 0212 

0 • 

2. 063 


-.3657 

.1579 

.0329 

1.1804 

.0307 

.3381 

1.630 


3884 

.1591 

.0896 

1. 1382 

.0359 

.7939 

1.013 


-.360 

.1257 

.1440 

.9504 

.0297 

1.224 

.630 


-.30 

.0907 

.1749 

.7858 

.0212 

1.493 

.404 


-.14 

.0343 

. 2063 

.5590 

.00774 

1.830 

.134 


0 

0 

,220 

.440 

0 

2.000 

0 


.50 

-.0897 

.2459 

.1793 

<-.0188 

2. 347 

-.257 


^ 2.00 

-. 2938 

,2829 

-.2933 

-.0586 

2. 837 

-.638 

- 0.8 

-0. 3088 

0.1215 

0 

1.0305 

0. 0172 

0 

2. 240 


-. 325 

.1304 

,0312 

1.0600 

.0210 

.3100 

1.828 


-.3285 

.1298 

,0436 

1,0499 

.0216 

.4194 

1.708 


-. 3285 

.1260 

.0681 

1.0185 

.0220 

,6245 

1.601 


-.325 

.1212 

.0827 

.9885 . 

.0216 

.7438 

1.396 


-.30 

.1017 

.1214 

.882 . 

• .0187 

1.058 

1. 138 


-.14 

.0355 

.1935 

.5781 

.00042 

1.712 

.692 


0 

0 

.220 

.44 

0 

2 . 0 (K) 

.619 


.50 

-.0837 

'.2678 

. 1676 

-.0138 

• 2.599 

. 199 


1.50 

-.2008 

.3179 

-. 1332 

-.0313 

3.263 

-.083 


2.00 

-.2522 

.3366 

-. 2517 

-.0388 

3. 502 

-. 166 

-0.4 

-0.246 . 

0. 0899 ’ 

0 

0.9087 

0. 00679 

0 

3. 041 


-. 2483 

.0894 

.0300 - 

.8968 

.00730 

.2941 

2. 679 


-.24 

.0826 

,0624r 

V . 8519 

.00696 

.5775 

2. 399 


-.20 

.0615 

.1210 

. . 7379 

.00554 

1.074 

2.034 


0 

0 

*.220 

.440 

0 

2.000 

1. 556 


.50 

-.0722 

.3019 

.1442 

-. 00573 

3. 042 

1. 185 


2. 00 

-. 1733 

.3924 

-.1713 

-.0118 

4. 628 

.759 

0 

-0. 1988 

0.0681 

0 

0.822 

0 

0 

4. 032 


-.16 

.0487 

.1051 

.7068 

" 1 

. 9480 

3. 094 


0 

0 

.220 

.440 


2. 000 

2. 691 


.50 

-.0602 

.3220 

.1232 

0 ' 

3. 436 

2. 298 


1.00 

-.8029 

,3556 

0 

0 

4. 317 

2.218 


1.60 

-.1002 

.3808 

-.0748 

0 

5 . 122 

2. 180 


2. 00 

-. 1064 

.3892 

1040 

0 

5. 565 4 

2.162 

1.0 

-0. 1295 

0. 0417 

0 

0.7280 

-0. 00803 

0 

6. 723 


-.10 

.0294 

. 09798 

,6476 

-.00644 

.8956 

5. 671 


0 

0 

.220 

.440 

0 

. 2.000 - 

5.187 


.30 

-.0334 

.3277 

,1558 

.00607 

3.602 

5. 493 


.50 

-.0735 

.3384 

.0755 

. 00588 

4. 315 

6. 012 


1.00 

-.0312 

.3065 

0 

.00338 

5. 644 

7. 850 


1.50 

-.0186 

.2382 

-.0114 

.00133 

6.662' 

11.125 


2.00 

-.0089 

.1663 

-.0088 

.00034 

7.627 

17. 105 


V 
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Now, if definitions (21) and (22) arc inserted, equation (26)' 
becomes 

A 

dX' 


-Ue ~ (^)=2[n(f/.+2)+/] (27) 


A similar procedure can be carried out with the energy 
equation, although it is not necessary for the calculation 
method herein presented. 

CORRELATION IN TERMS OF n 

If, in the definition of (eq. (IG)), S is set equal to zero, 
the form factor then becomes the same as that of Thwaites 
and equation (26) becomes Thwaites’ momentum equation. 
If, in this case, physically valid relations Hir(n) and l(n) 
can be established, the equation can be integrated and the 
problem is solved. The assumption of the form of the veloc- 
ity profile serves this purpose, and the resultant procedure is 
the well-known Pohlhauscn technique. (Kalikhman (ref. 
16) was the first to carry the same approach over to the case 
of the thermal profile.) Thwaites used the more direct ap- 
proach’ by determining whether universal relations Htr(n) 
and l{n) could possibly be established from the well-known 
exact solutions of the boundary-layer problem. An exami- 
nation of these solutions proved that for favorable pressure, 
gradients a single relation for each of these quantities could 
be' established with a fair degree of accuracy, but for adverse 
pressure gradients the relations departed from each other 
considerably, as indicated for l(n) in figure 1 (taken from 


The concept of correlation is herein extended by the 
following major assumption: Fo r the compressible lami nar 
boundary Jayi'T wilb lu'nt, I.i-n.nsb'i- nfi’nsa tb(\ siiffncp^ the 
skim=fr4clio»--4wd-Jm4rL44'-ausf(a:— paxa jnctcrs I and r can be 
correlated only in terms of the parameters 7i nnd Sw, derived 
l^’ ifrom the exact solutions of reference 1^ 7~ - 

— -I ■ It is thus implied tlial tlic solutions of reference I'J 

adequately represent tbe general boundary layer, although 
-thov- were derived for Falkncr-Skan type f low. * 

As in all firsf^-order boundary-fayer theories, the pressure 
distribution (and consequently the external velocity dis- 
tribution) is assumed to be known. Then the utility of the 
correlation may be stated as follows: If n i s known at a 
given point on the surface, 0<r (and hence d) canTmmcdiatcly 
be obtained from equation (22). If l{n) is a known function 
for the specified wall temperature, the wall shear is im- 
mediately obtainable from equation (21). Similarly, if 
r{n) is known, the heat transfer can bo found from equation 
(23). 

If the postulate of correlation in terms of n and Sw is 
admitted, equation (27) becomes 


m.z 


.--Hortree (ref. 26) 
.-Howorth (ref 25) 
-Thwaites (ref. 21 ) 
-Pohihousen (method 
of ref. 15) 



~u. 


n 


where 


A/ 

dx 


> 


■-N(n,Sw) 


N(n,Sw)^2[n{H„-i:2)+l] 


(28) 


(29) 


Figure 1. — Slsin-friction correlation of Thwaites (ref. 21) for incom- 
pressible flow over an insulated surface. 


This is the fundamental equation of the present method. 
Its solution, resulting in a determination of njxj , is the first 
stage in solving for the boundary-layef~%Iiaracteristics. 
Then the function l{v,Sw) is used to determine the wall 
shear, and the funcfioTrr'(7i-,N„) is used to determine the 
heat transfer. 

EVALUATION OF CORRELATION FARAMETERS 

, . ' The quantities I, n, and r definqd in, cquations (21) to 
;■ ?i23), as evaluated from the solutions 'of reference 12, .arc 
] i| listed in table II. An alternate parameter to r for the de- 
termination of the heat transfer is the Reynolds analogy 


ref. 21). Since the correlation technique requires the selec- 
tion of a single relation between all boundary-layer quanti- 
ties regardless of the history of the development of the spe- 
cific boundary layer under 9 onsideration, the assumption of 
a single relation for all boundary layers is not exactly valid. 
Thwaites chose a relation for l(n) that would match the 
Howarth solution (ref. 25) at separation. The Pohlhauscn 
solution predicts separation to occur much later than any 
of the other solutions in the sketch, while the Hartree solu- 
tion (ref. 26) predicts separation to occur earlier than the 
Howarth solution. Stewartson (ref. 24) has indicated that 
under certain emulitions Howarth’s solution would jircdict 
separation too li)^ 

"With the exception of the described differences in the 
relation between boundary-layer quantities (due to tlic • 
selection of the solutions of ref. 12 for the evaluation of the 
correlated quantities), the method of correlation to be 
presented herein will contain as special cases the method of 
Thwaites for incompressible flow and the results of Rott 
and Crabtree for compressible flow over insulated surfaces. 


parameter, defined as relating the heat transfer to 


the skin friction. Because this (arbitrarily chosen) parameter 
is herein determined from solutions for a Prandtl number of 

I. 0, it will be denoted as, it appears, in table 

II. It is related to r by 

2Sanl 


/ Cfllew\ 
\ Nu 

The parameters I and 


(30) 


as ' functions of n 


and Sw are plotted against n in figures 2 and 3, respectively. 
The solid portions of the curves represent the solutions of 
reference 12. The reversal of the curves for »S'„=1.0 is 
associated with the velocity overshoot phenomenon discussed 
in reference 12. 




Dimensionless shear parameter. 
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METHOD OF SOLUTION 

. The solution of equation (28) obviously depends on a 
knowledge of the term N(n,Su,). This quantity was evalu- 
ated from equation (29) and the associated formulas for 
l,n,andHtr. Th e results are shown in figure 4. It is to be 
noted that the curve for 0 is nearly a straight? line, and 
for all negative values of Su, the curves depart only slightly 
from this condition except in the range near separation 
where the curves become double-valued. For Sa'^0 (hot 
wall) the curve is essentially straight except in the region of 
strong favorable gradients. 

The examination of these curves and of equation (28) 
suggests two methods of determining the correlation number 
n. The first method is that of solving equation (28) by 
numerical integration. The necessary numerical procedure , 
however, is tedious, since it involves integration of a firs t- 
order nonlinear nonhomogeneous ordinary differential co ua- 
tion. Because a simpler method is available when the 
surface is isothermal, no numerical integration procedure 
' will be presented here in detail. However, some integration 
relations are stated in appendix B. 

The second method, applicable when the surface tempera- 
^turo is constant (or, presumably, nearly constant), will bo 
I termed the “linear method.” This method uses the nearly 
linear shape of the curves of N against n for constant 8a. 
lilt directly corresponds to the procedure of Thwaites for 
' Uncompressible flow and to that of Kott and Crabtree for 
'compressible flow over insulated surfaces. - The curve of 
N against n for a given is assumed represented by 



Correlation number, n 


Figure 4. — Correlation of dimensionless momentum equation. 


N=A+Bn (31) 

If equation (31) is inserted in equation (28), a simple 
linear first-order equation results, which has for its solution 


n 


= -AUr^U. 


■ 4 : 


U?-^dX 


(32) 


If equation (32) is transformed to physical quantities by 
using Stewartson’s transformation, there results for two- 
dimensional flow 




(33) 


P' i!: 

t 


where K= 


3y — 1 


■ 2 ( 7 - 1 ) 

less pressure gradient P' is given by 
j dp^ i 

p,_ ' dx _p c dx 


L is any fixed length, and the dimension- 


L du. 


peUl yMl dx 


(34) 


/ 

[ The left member of equation (33) has been arranged in a 
[i' form convenient for later use. 

The determination of the coefficients A and B is as follows: 
If the straight line (31) is chosen to pass through the correct 
value of N at zero pressure gradient (n=0), then H.=0.44 
I independent of S'„. In this case, only B is affected by the 
! presence of heat transfer. Figure 5 shows the values of 



Figure 5. — Variation of B for use in linear method of determining 
correlation number. 
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B{Sm) for the following choices of matching conditions: 

( 1 ) The line (31) goes through the point corresponding to 
N=0 (y.e=0) in order to match the conditions at a stagnation 
point. Thus, for two-dimensional flow (iV=0), 

— 0.44/risp. 

( 2 ) The lino (31) coincides with the tangent to the N(n) 
lines at n =0 (stnall pressure gradients). 

(3) The line (31) is selected to give good over-all agree- 
ment for unfavorable pressure gradients. 

If a better match with the curves of figure 4 is desired in 
calculating n for certain ranges of pressure gradient, a 
tangent line to the curve of N against n may be chosen at a 
desired' value of n. For instance, in considering the flow in 
(the vicinity of a two-dimensional stagnation point with 
<S'k=— 0.8, the tangent line through iV=0 has the coefficients 
4=0.372, B=2.53. The value of B is quite sensitive to 
jhe matching assumption, especially in the region iSa^O 
\ (fig. 5). However, the final value of n(x) is somewhat in- 
, sensitive to the value of B, since 4hc-terms involving B in 
equation (33) appear both inside and outside the integral in 
a compensating manner. The accuracy of the method de- 
creases, of course, in regions where the plots of figure 4 have 
large curvature. 

The calculation procedure is as follows: '.Values of A and B 
are chosen for use in equation (33) either from figure 5 or 
from tangent-line considerations. The integration is then 
performed by using a suitable integration rule and a proper 
step size. It is recommended that the step' size chosen be 
as small as practicable in order to obtain results which are 
reasonably smooth. In some cases (e. g., near the boundary- 
layer origin) it may bo advisable to perforih the integration 
by obtaining Taylor series expansions of the integrand in 
the variable {xjL). Then the integration can bo carried 
out in closed form, corresponding effectively to zero step 
size. 

There are two possible starting conditions in a boundary- 
layer calculation: (a) that of a sharp edge, that is, 0 = 0 , 
n=0; or (b) the stagnation point, where Ue—Ue=0. In 
using the linear method, the starting conditions are auto- 
matically satisfied if the chosen lino (31) goes through the 
correct starting point N{n). Thus, if matching condition 
(a) is used, both possible starting conditions can be satisfied, 
since the corresponding line (31) passes through, the curve 
from the exact solutions at both n=0 (V=0.44) and, for 
the two-dimensional flow, at A ^=0 («,,„=— 0.44/H). Values 
of n for stagnation-point flow taken from figure 4 are shown 
in figure 6 . 

The corresponding relations and procedure for axially 
symmetric flow based on Mangler’s transformation (ref. 27) 
are presented in appendix C. 

It is sometimes helpful to have an analytical expression 
for the initial variation of correlation number as a check on 
the numerical calculations. The initial variation of n with 
X for the various starting conditions, as represented by the 

derivative ! is discussed in appendix D. 



Figure 6. — Effect of wall temperature on correlation number . for 
^ D P tk* various pressure gradients. 

BOUNDARY-LAYER CHARACTERISTICS 




Once the correlation number n is determined as a function 

I — -■ "" / (j \ 

; of X, it is possible to obtain I and from figures 

. 2 and 3, respectively. Then, the local skin-friction coefficient 
and heat transfer are easily calculated from the following 
relations, which apply to both two-dimensional and axially 
symmetric flows: 

, Re^=^, Nu= 


>!4 


From the definitions 


2 Pw^e 


iav 


and from equations ( 21 ) and ( 22 ), it follows that 





(3f.) 


It may be noted that, at a stagnation point, equation (35) 
reduces to 

- 7 = / (3G) 

“V 'P'sp 


Cf-^Rsu 


Once Cf is determined, the heat transfer may be calculated 
from curves of the Reynolds analogy parameter against corre- 
lation number of figure 3, by using the relation 


Nu Cf^Reg 

^W~(CfR^ 

■ V Nu Jpr^l 


(37' 
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In utilizing equations (35) to (37), it is useful to have the 
initial values of the parameters. These values are listed in 
table I. 

The calculations thus far have been for Pr= 1,0. The 
effect of Prandtl num.ber on skin friction is small and is there- 
fore usually neglected. It can be seen from table I that, for 
stagnation-point flow, the maxim.um. difference in the quan- 
tity between solutions for'Pr=1.0 and Pr=0.7 is 

about 7 percent. With regard to heat transfer, Tifford and 
Chu (ref. 28) have found, from solutions with constant fluid 
properties, that the effect of Prandtl number on heat transfer 

can be accounted for by multiplying f ) by (Pr)“. 

Using this approxim.ation, equation (37) becomes 



Values of a suggested in reference 28 arc as follows: For 
small pressure gradients, a=}i; for large adverse pressure 
gradients, a=%; and for extrem.ely favorable gradients, a=K. 
Squire (ref. 29) has indicated that a=0.4 is adequate for 
stagnation-point flows. Kecently obtained solutions (ref. 30) 
of equations (16) of reference 12 for |3=1, Pr=0.7, and 
“ show that this typo of correction m.ay bo adequate 
for all compressible boundary-layer calculations. 

[ It should be noted that, in the definition of Nusselt' 
! number, the temperature difference in the denominator was 
j assumed to be (taw—tw)- Since for Pr=l the reeovery 
j temperature is U, the present calculations (based on those of 
ref. 12) can give no indication of the adiabatic wall tem.pera- 
ture for Pr?^\. For a first approximation, it may be reason- 
able to calculate taw by using a temperature recovery factor 
of (Pr)'^^. This is the well-known expression for recovery 
factor for the case of high-speed flow with zero pressure 
I gradient. The adequacy of its application to high-speed 
i flows with largo pressure gradients is not well established. 

The physical momentum thickness is determined from 



The displacement thickness 5* may be calculated by using 
, the following simple expression for the ratio of displacement 
■ , thickness to momentum thickness: 


I Ml{H, r+1) (40) 

In reference 23, this expression was derived for flows over 
insulated surfaces with Pr=l. Equation (40) is valid for 
noninsulated surfaces as well. The dependence of Htr upon 
wall tem.perature and ‘n is presented in figure 7. With large 
amounts of cooling in favorable-pressure-gradient flows, it is 
seen that negative form factors (and hence negative dis- 
placement thicknesses) result. This occurs because the sur- 

I face cooling produces an increase in density near the wall, so 

i 

J 

I • ■ ' ( 

I ■ f 

( 

; 



Correfation number, n 
Figure 7. — Form-factor correlation. 


that there is more m.ass flow per unit flow area wdthin the 
boundary layer than in the external flow. 

The over-all thickness S of a boundary layer calculated 
from, exact solutions is a quantity not uniquely defined. Its 

value depends on the value of the velocity ratio — that is 

chosen to represent the outer edge of the boundary layer. 
However, for a given value of the velocity ratio (less than 1) 
there is a single value for thickness. The ratio of this over-ail 
thickness to the momentum thickness is given by the expres- 
sion 

(41) 

The quantity ^ is that for low-speed flow (Me=0) which 

^tr 

has been evaluated from, the solutions of reference 12 and is 
presented in figure 8 as a function of correlation number n 

for ^=0.995. 

Ua • . 

The mass flow in the boundary layer is related to the differ- 
ence (5-5*) and may be obtained by subtracting equation 
(40) from, equation (41), which finally results in 
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The quantity is 5.35 for the flat plate and, as can 

be seen from the information in figures 7 and 8, varies from 
about 4.5 for adverse pressure gradients to about 6.5 for 
favorable pressure gradients. ^ 

LIMITS DUE TO AVAILABLE EXACT SOLUTIONS 

The curves of figures 2 and 3 are terminated at the boundary- 
layer separation point on one end and at the value of the 
correlation number at the other end, which corresponds to 
the most favorable pressure gradient that can be represented 
by the Falkner-Skan type solutions of reference 12. How- 
ever, one may conceive a flow that contains values of n more 
negative than the latter limit. For example, if a boundary 
layer is allowed to grow on a flat plate until the momentum 
thickness has an arbitrary value and is then subjected to a 
favorable pressure gradient («6j,>0), equation (22) shows 
that n might be made arbitrarily negative. Similarly, large 
negative values of n may arise if a strong favorable pressure 
gradient is maintained for some distanee along a surface, 
such as in a nozzle. Under such circumstances it would be 
necessary to extend the curves of figures 2, 3, and 7 in order 
to use the correlation m.ethod presented herein. Sueh exten- 
sion is not possible by m.eans of Falkner-Skan type solutions. 

It is thus apparent that there are types of problems for 
which the presently available correlations are not adequate. 
To establish correlations in that regime would require exact 
boundary-layer solutions for large favorable pressure gra- 
dients of a form different from that of reference 12. Such 
solutions are not known to the authors at this time. The 


extended portions of the curves in figures 2, 3, and 7 for 
S„——0.8 were drawn to enable the completion of example 
II in the following section. The adequacy of this extension 
is not known. 

EXAMPLES 

An important test of the method developed is the com- 
parison of the final results for practical problems with the 
findings of other theories or with experimental results. 

I. SUPERSONIC SURFACES 

The linear method for determining the correlation number 
n is applied to the calculation of skin friction and heat 
transfer for the two supersonic surfaces at Mach number 3.0 
calculated in reference 6. These surfaces arc shown in 
figure 9. 



x/L 



A comparison is made, in the following table, between the 
results obtained by using the linear method and those ob- 
tained by using Low’s perturbation method with Pr=0.72 
(ref. 6) . The comparison is made for a hot wall, an insulated 
wall, and a cold wall at xjL—l (see fig. 9). A value of a of 
1/3 was used in these calculations. 


^^radlent 

Heat-transfer 

condition 

{Cf-^ReuJx/L^i 

• {Nul‘^ReJ)g/L^i 

Linear 

method, 

Pr=l 

Low 
(ref. 6) , 
Pr=0.72 

Linear 
method, 
corrected 
for Pt 

Low 
(ref. 6), 
Pr*=0.72 

Favorable. 

Cold wall (S,.=-0.9). 

Insulated wall 

Hot wall (S«>s=0.61)... 

0. 680 
.877 
1.031 

0. 679 
.876 
1.048 

0.270 

,284 

0. 271 
.318 

Adverse... 

Cold wall (S«, = -0.9)_ 

Insulated wall-... 

Hot wall (S„=0.43)... 

0.645 

.406 

.256 

0. 661 
.475 
.386 ■ 

0.320 

.274 

^ 0.307 
.277 


, A comparison of values indicates agreement of skin friction 
within 2 percent in the case of a favorable pressure gradient. 
For the adverse-pressure-gradient cases, reasonable agreement 
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is obtained for the cold wall, although for the insulated and 
hot walls a largo difference is obtained. This difference is 
due essentially to the fact that, in the case of an adverse 
pressure gradient, the solutions of Low (ref. 0), which re- 
semble the series-typo solution of Howarth (ref. 25), depart 
from Falknor-Skan typo solutions such as that of Hartrcc 
(ref. 2G) (o. g., fig. 1). An important consideration for the 
case of the heated surface with an adverse pressure gradient 
is that the flow is closer to separation than aj)j)cars iiermissi- 
blc for a theory based on small pressure gradients such as 
that of reference G; therefore, for this case the present calcu- 
lation may bo the more reliable. Some of the diirorence in 
the preceding table may also bo a Prandtl number effect. 

Good agreement is also obtained for heat transfer except 
for the case of the heated surface with favorable pressure 
gradient. Some of that difference might also bo a Prandtl 
number effect. 

II. AXISYMMETRIC CONVERGENT-DIVERGENT ROCKET NOZZLE 

The second example, that of a rocket nozzle, is one involv- 
ing both largo pressure gradients and heat transfer. The 
nozzle chosen is illustrated in figure 10. It has a 25° half- 
angle convergent section and a 15° half-angle divergent 
section. The combustion-chamber stagnation pressure is 
assumed to bo 500 pounds per square inch absqlutc, the 
stagnation temperature is taken as 4000° K, and the Prandtl 
number is assumed to bo 0.78. The nozzle wall is assumed 
cooled to a uniform temperature of 800° R, which corre- 
sponds to (S'„=— 0.8. For the assumed 3-inch throat di- 


ameter, the rocket has a nominal thrust of 5550 pounds for 
7=1.3. Local static conditions along the nozzle wall were 
obtained using one-dimensional flow relations. 

The calculation was jierformed by the linear method with 
Al=0.372 and R=2.53. In order to eliminate the effect 
of step size in the initial portion of the integration, a series 
expansion of the integrand was used for 6 < (a;/Zl) <(0.5. 
For (x/Zl)(>0.5, the stei) size taken was 0.1. The resulting 
variation of n in the nozzle is also shown iniigurc 10. It is 
seen (hat, in a iiortion of the nozzle including the throat 
(1.2<((x//>^)<C3.5), values of n'are obtained which require 
use of an o.xtoiulcd correlation curve as discussed earlier, in 
order to calculate skin friction, beat transfer, and displace- 
ment thickness. No c.xtrapohition is needed to obtain 
momentum thickness, since the momentum thickness is 
related to n through equation (39) . 

The calculated local heat-transfer rates as well as displace- 
ment and momentum thickness arc shown in figure 11. It 
is seen that largo rates of cooling are required in the neighbor- 


Nozzle throat 


/ \ 
/ 1 \ 




o Exxt solution for 

stagnation-point flow . 

Exfended correlations 

used 
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hood of the stagnation point and the nozzle throat. If the 
cooling were to become insufficient, these seem to be the most 
likely locations of failure. The required local cooling rate 
(o maintain constant wall temperature decreases sharply 
beyond the throat of the nozzle. 

In the absence of more appropriate information, it has 
been customary in recent years to use flat-plate hoat-transfer 
relations in estimating heat transfer in a nozzle. The use 
of such a relation in the current problem (A%/VA’e„=0.,‘{05 
for Pr=0.78) yields values indicated in figure 10. It is 
seen that for this problem the flat-plate relation seriously 
underestimates the amount of cooling required over a largo 
part of the nozzle. This illustrates the importance of con- 
sidering the effects of pressure gradient on heat transfer. 

The momentum thickness is seen to reach a minimum 
value at the nozzle throat. "Idio displacement thickness 
is a small positive quantity at the stagnation point but is 
negative for most of the convergent section as well as in 
the vicinity of the throat.^ The over-all thickness (not 
shown on the figure) varies from 6.1 to 7.7 times the momentum 
thickness in going from the initial stagnation point to the 
nozzle exit. 

The use of different values of A and B in performing the 
numerical integration would have approximately the fol- 
lowing effect: With the values. .4=0.44 and .5=3.0 (from 
fig. 5), the momentum thickness would be about 10 percent 
smaller in the vicinity of the throat than the values in figure 1 1 
and would be within 5 percent of the presented values over 
the rest of the nozzle. With 4.=0.335 and 5=2.34, the 
momentum thickness at the throat would be about 6 percent 
larger than the value presented. The effects of varying 


* This unusual result produces the interesting possibility that, for a rocket nozzle with 
cooled walls and viscid flow, a mass discharge coefficient based on throat area, generally 
assumed to be less than 1 because of boundary-layer “blockage" at tne throat, may actually 
exceed 1. A distinction exists between this phenomenon and that of negative momentum 
thickness (refs. 1, 10, and 12) associated with velocity overshoot. 


A and 5 on skin friction and heat transfer would be less 
than .3 percent at the throat. 

CONCLUDING REMARKS 

The ap])lication of Stewartson’s transformation to the 
compressible laminar-boundary-layer equations with heat 
transfer yielded a simple first-order system of ordinary 
difl’erential ccpiatiojis, the first of which is very similar to 
the Karman momentum integral. Dimensionless shear and 
heat-transfer parameters were defined. The assumption of 
correlation of these parameters in terms of a momentum 
parameter resulted in a complete system of relations for 
calculating skin friction and heat transfer. Knowledge of 
velocity or temperature profiles is not necessary in using 
this calculation method. Procedures for the calculation 
of the longitudinal distribution of correlation number arc 
presented, which include as special cases the method of 
Thwaites and that of Rott and Crabtree. The dimensionless 
parameters introduced herein were evaluated from the exact 
solutions of reference 12. 

Calculations of an example involving small pressure 
gradients have shown reasonable agreement between this 
method and the perturbation method of reference 6 over 
the same range of Mach number, pressure gradient, and ' 
heat transfer. 

The method is also applied to the calculation of heat 
transfer and displacement thickness in a highly cooled, 
convergent-divergent, axially symmetric rocket nozzle. The 
results of this calculation show that high rates of heat transfer 
are obtained at the initial stagnation point and at the throat 
of the nozzle. Also indicated are negative displacement 
thicknesses in the convergent portion of the nozzle; these 
occur because of the high density within the lower portions 
of the cooled boundary layer. 

Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, February 1, W55 
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APPENDIX A 


SYMBOLS 


A constant from 

a sonic velocity ; 

B constant from N=A-\-Bn 


Cf local skin-friction coefficient, Gj= 

specific heat at constant pressure 
nozzlo-tliroat diameter 


d ■ 

IS 

H 

Htr 

h 

K 

k 

k 

L 

I 

M 

m 

N 


n 


P 

II 


S 


2t„ 

Pw^e 


fA U 

convection tliickness, E= Sjj-dY 

Jo Ue 

form factor, H=8*jd 

S* 

physical form factor for Me~0, IItr=^ 

“i' 

enthalpy 

37-1 

2(7-1) 

thermal conductivity 
Sutherland’s constant 
arbitrary length 

dimensionless shear parameter, 

Mach number 

exponent from Falkner-Skan external velocity distri- 


bution Ue=CX” 
momentum parameter, iV=2[n(//„-l-2)-t-Z] 


Nu Nusselt number, 


ta 


correlation number, n= 


> Zw 

■ Vo 


L dp, 

P' dimensionless pressure gradient, P'= 


Pr Prandtl number, Pr- 


JlCj, 


static pressure 

radius of axisymmctric body 

pvMe'X 


2Ze„ Ecynolds number, Ilea=- 


Pu 


heat-transfer parameter, 


enthalpy fimction, S=^—l 
ho 


t static tcmi)craturc 

taw adiabatic wall temperature 


U 


transformed longitudinal velocity, U=u—=^y 


u longitudinal velocity component 

V transformed norihal velocity, V=—^x 

V normal velocity component 

T • Ct 

X transformed coordinate along surface, -X’= \ — — dx 

Jo Poao 

X coordinate along surface 

V transformed normal coordinate, Y=— \ —dy 

aojo po 

y normal coordinate 

a exponent of Prandtl number in Keynolds analogy 

parameter 

2m 


|8 pressure-gradient parameter, /3= 

7 ratio of specific heats 

5 over-all thickness 

8* displacement tliickness 

6 momentum thickness 

■y ( m/ Mo ) / Zo~t~ ksu \ l ^w 

(t/to) \tw-\'ksu/\ to 
p dynamic viscosity 

V kinematic viscosity, v=pjp 

p mass density 

bu 


m-\-l 


shear stress, r = p 


by 


\f/ stream function, eq. (6) 

Subscripts: 

■ e local flow outside boundary layer (external) 

s local stagnation value 

sp stagnation point 

tr associated transformed quantity 

w wall or surface value 

0 frcc-strcam stagnation value 

1 initial value ^ 

A coordinate used as subscript denotes differentiation 

with respect to the coordinate. 
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APPENDIX B 


NUMERICAL INTEGRATION METHOD 


METHOD 

The most direct method of solving equation (28) is by 
numerical integration, using the calculated curves of iV(R,5'u,) 
for determination of the righ(i member. 

An integration procedure may bo simply indicated by 
direct integration of equation (28). The resultant integral 
equation can be written: 

For two-dimensional flow, 




For an axially symmetric closed body, through Mangler’s 
transformation (ref. 2G), 


t/, 


n=~ 


X 

IP 


NIP 

Jo u. 


dX 


(B2a) 


Since the integrands contain N{n), which is unknown at 
this point, no simple evaluation is possible. In fact, these 
equations are actually only a condensed notation of the pro- 
cedure to be followed. The integration must be carried out 
piecemeal, alternating with determination of the left member 
of the equation and iterating for accuracy. 

The necessity of working with the transformed coordinates 
can be eliminated by considering the Stewartson trans- 
formation from Ue to We- For example, in physical coordi- 
nates equation (B2a) becomes 



where L is any fixed length, and the dimensionless pressure 
gradient P' is given by 

jr dPe L dPe 

p, dx Pe dx _ (B3) 

~ PeW2 ~ 7M| 

Similarly, if the iscn tropic relation p/p'’' = constant is used 


in equations (31) and (32) and if the Stewartson transfor- 
mation is applied, there results: 

For two-dimensional flow. 



For an axially symmetric flow (closed body). 



where K-- 


37-1 

2(7-1)’ 


INITIAL VALUES 


(B4) 


(B5) 


When the numerical integration m,cthod is used, certain 
considerations are necessary in order to start the solution 
properly. There are two possible starting conditions: 

(1) sharp edge or pointed body, where 6=0 and n = 0, and 

(2) stagnation point, where Ue=Ue=0. 

In the case of a boundary layer starting from a stagnation 
point, the initial value Wi of n is determined from the condi- 
tion Ue=0, Z7„=constant in equation (28). For two- 
dimensional stagnation -point flow, the Hartreo pressure- 
gradient parameter ^ is equal to 1.0. Since, for the Falkner- 


Skan type flow considered, N= 


2(|3-1) 

/3 


it is seen that 


N=0 at a two-dimensional stagnation point. This fixes n 
at the values shown in figure 6, which were obtained from 
figure 4. For axisymmetric stagnation-point flow over a 
closed body, ff=}i (ref. 31), so that Wi = — 2wi. The values of 
ni for axially symmetric stagnation-point flow over a closed 
body as obtained from figure 4 arc also shown in figure G. 
For the stagnation-point flow over the blunt lip of an open 
axisymmetric body, risp can bo shown to have the two- 
dimensional value. 
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APPENDIX G 


LINEAR METHOD FOR AXISYMMETRIC FLOW 


For axisymmetric flows, the following equation, which is 
equivalent to equation (33), is obtained by application of the 
transformation of Mangier (ref. 27): 


^(Fo) "7^ Jo (f) Mp (i) (Cl) 


P' — 
te 


where l?=R{x) is the radius of the body at station x, 
K =^ — and 


2(7-1) 


P'= 


T djpe 
dx 


yPeMl 

In cvaliiating the coeflicicnts A and B of the straight line 

(31) 


A may ho chosen as 0.44 so that the line (31) passes through 
the correct value of A” at n=0. The choice of B may bo 
made so that the lino (31) goes through the point where 
iV= — 2n in order to match the conditions at an axisymmetric 

/0.44 \ 

stagnation point. In that event, B=—{ [-2 )• 

\ Xl:fp / 


For 


achieving better agreement with the curves of figure 3 in 

certain ranges of pressure gradient, a tangent line to the’ 
curve of N against n may bo chosen as was indicated for 
two-dimensional flow. 

The three possible starting conditions arc: (1) For a 
pointed body, 0=0, n=0. (2) For a stagnation point on a 

closed body, Ue—Ue=0 so that N——2n. Values of n,p for 
this axisymmetric stagnation point are shown in figure 5 and 
indicated in table IT for |3=H. (3) For a stagnation point 

on the blunt lip of an open axisymmetric body, it can be 
shown using the axisymmetric form of equation (28) that 
N=0 so that Wsp is that for two-dimensional flow. 


APPENDIX D 

INITIAL VARIATION OF CORRELATION NUMBER 


It is sometimes helpful to have an analytical ejiprcssion for 
the initial variation of correlation number as a check on 
the numerical calculations. The following expressions for 

are determined from equation (28) : 

TWO-DIMENSIONAL FLOW 

Sharp edge: 

(p) =_o.44(^) 

\dX/ a: = o \ He /z = 0 \ ie/ 

Stagnation point (blunt body) : 


(Dl) 


/ ^ risp I 

\dz/j = 0 1 I \^‘z)s 


(D2) 


AXISYMMETKIC FLOW 

For axisymmetric flow the initial derivatives must be 
evaluated from the following equation : 


Ue d /nR‘-\_ 
IP dX (f/e 


>) 


(D3) 


Closed body. — 
Pointed nose: 


, (p) =-0.147 (^) 

\(iX/x = o \ '^e /x = Q \tc/x=0 

Stagnation point (blunt nose) : 


Open body.— 
Sharp lip : 


.(P\ =_o.44(^) (p) 

\(ix / x = p \ / X=0 *=0 

Stagnation point (blunt lip): 

\dxJi:=o X I \P/spJ 


(D4) 


\ dlT Jsp 


(D6) 


(D7) 


/ 
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